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Abstract 

The effective QCD-interaction as obtained from the front form Hamiltonian by DLCQ is 
Fourier transformed on the Retarded Schrodinger Equation for to describe the constituents of 
physical mesons. The crucial point is the use of a running coupling constant a s (Q 2 ), in a manner 
similar but not equal to the one of Richardson in the equal usual-time quantization. Fixing the 
running coupling constant at the Z mass, the only parameters are the flavor masses. Without 
the top quark one needs thus 5 parameters to calculate the physical masses of 30 pseudoscalar 
and vector mesons, consistently within the same model. Applying variational methods to a 
caricature of the model, the biggest technical challenge is the solution of a cubic algebraic 
equation. - In view of an oversimplified model and a very simple technology, the agreement 
with the empirical data is much too good. 
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1 Introduction and Motivation 



One of the most outstanding tasks in strong interaction physics is to calculate the spectrum and the 
wavefunctions of physical hadrons from Quantum Chromodynamics. The method of 'Discretized 
Light-Cone Quantization' (DLCQ) [0 has precisely this goal and has three major aspects: 

(1) a rejuvenation of the Hamiltonian approach, 

(2) a denumerable Hilbert space of plane waves, 

(3) Dirac's front form of Hamiltonian Dynamics. 

In the front form, using Dirac's terminology |Q, one quantizes at equal 'light-cone time' x + = 
x° + x 3 , as opposed to the instant form where one quantizes at equal 'usual time' t = x°. The front 
form carries many names like for example light-cone quantization Q, see also Q. The approach 
has some unique features as reviewed in Q: The vacuum is simple, or at least simplier than in the 
instant form, which in turn implies that the groundstate of the free theory is also the groundstate 
of the interacting theory. The relativistic wavefunctions transform trivially under certain boost 
operations [Q, |4|. Both are in stark contrast to the conventional instant form. Discretization, 
that is the use of periodic or antiperiodic boundary conditions, provides additional length scales 
as useful regularization parameters in a formalism like a quantum field theory which otherwise is 
littered with divergencies ab ovo. In fact discretization (or working 'on a circle') is more than only 
a useful tool for calculations, rather it provides a systematic scheme to formulate the theory at all. 

Over the years, the light-cone community ||] has made much progress. Calculations have been 
done Q to verify agreement with other methods particularly lattice gauge theory. We have learned 
that zero modes should not be ignored ||, since they can be important carriers of quantum structure. 
Other calculations in 1-space and 1-time dimension like the dimensionally reduced models [Q, U 
have an interest of their own. We have learned that there are physical and non-physical gauges, 
even something like a wavefunction of the vacuum ^ . Renormalization is a difficult problem within 
a Hamiltonian theory and worth the effort of the best || [T^] . 

What was not shown thus far, however, is why all that is useful. Particularly one lacks the 
contact to phenomenology beyond the perturbative regime We believe that one needs more 
work like or like Jl^, [l3j which relate formalism and actual experimental numbers. It should be 
possible to combine rigor with simplicity! The present work will hopefully be a useful contribution 
to the extent that we shall derive a rather simple, conventional second order differential equation for 
the wavefunction in configuration space. Amazingly enough, it seems possible to describe quarks 
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in hadrons by a non-relativistic Schrodinger equation to account for the correct retardation and 
still being in the front form of Hamiltonian Dynamics applied to Quantum Chromodynamics. 

Solving this 'Retarded Schrodinger Equation', or better a caricature of this equation, by para- 
metric variation, we end up with the masses of all pseudoscalar and vector mesons. In comparison 
with the experimental data |l4j| , they are no worse than conventional phenomenological models 
|15, [Tql , or predictions from heavy quark symmetry [17|, or even predictions based on lattice gauge 
calculations fl8| , O , 20]. In view of the simplicity of the methods applied we regard this as a success 



of the front-form approach, particularly since we can and will do better in the future 31]. 



2 The Model: The Effective QCD-Hamiltonian from DLCQ 

In Discretized Light-Cone Quantization one addresses oneself to solve an eigenvalue problem 

H LC \ipi) = EM) . (1) 

The operator Plc = P^Pfi is a Lorentz-scalar with the dimension of a <mass> 2 . The eigenvalues 
have the same dimensions. We therefore set Ei = Mf and interprete Mj as the 'mass of the state 
i\ In DLCQ one works in momentum representation, and thus three of the four components of 
P^ are diagonal, namely the three total (light-cone) momentum operators P + and P±. The fourth 
component, the total (light-cone) energy operator P~ is very complicated and off-diagonal. The 
four components of energy-momentum P^ mutually commute; their matrix elements can be found 
explicitly in Q. The eigenvalues of P + and P± are 

P+ = £>/ , and Pl = £4, , (2) 

j j 

where j runs over all partons (particles) in a Fock state; each parton has a four-momentum denoted 
by kj = (kj,kj,k±j) and sits on its mass shell (k^k^)j = mj. Based on the boost-properties of 
light-cone operators Q one always can transform to an intrinsic frame, where Pl = and thus 
P^P^ = P + P~ . Since P + is diagonal from the outset, the diagonalization of P~ and of Plc 
amounts to the same problem. 

The Hilbert space for diagonalizing P~ is spanned by all possible Fock states which have fixed 
and given eigenvalues of P + and Pj_ = (0, 0). They can be arranged into 'sectors' n = 1, 2, . . . , N, 
corresponding to their parton composition. We convene to call the sector with one quark and one 
antiquark as the sector with n = 1. The sector with one quark (q), one antiquark (q), and one 
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gluon (g) shall be denoted by n = 2, and so on. It is peculiar to DLCQ that the number of sectors 
is limited for any fixed value of P + , or of the 'harmonic resolution' K = irP + /L [0, since the 
longitudinal momenta kj~ are all positive and non-vanishing numbers. Both K and N can thus be 
arbitrarily large but are strictly finite. The eigenvalue equation, Eq.(||), can thus be written as a 
matrix equation for a finite number of block-matrices {n\H\n') 7 i.e. 



( 



<l|tf|l) 
<2|#|1> 



(1\H\2) 
(2\H\2) 



(1\H\N) 
(2\H\N) 



\ 



( 



(1| 

(2| 



M 2 



(1| 
(2| 



(3) 



^ (N\H\1) (N\H\2) ... (N\H\N) J \ (Nty) J \ W) J 

The projected wavefunctions are denoted by (n\ip). The probability amplitude for finding only a 
qq-p&ir with no gluons, in a pion for example, is denoted by (l\ip) = V^/ttj an d correspondingly 
for the higher sectors. The Hamiltonian matrix Eq.(|3|) has a finite number of sectors, but thus far 
each sector has an infinite dimension, since the transversal momenta of each parton can take any 
positive or negative value, including zero. The number of states is however finite if one admits only 
those transversal momenta whose absolute values are smaller than a 'cut-off'. The regularization 
of Brodsky and Lepage, i.e. 



X 



with Xn 



(4) 



does that in a Lorentz-invariant way. The cut-off mass A n can but need not depend on the sector 
number. In the qq sector, we shall use below A 2 = A 2 + (mi + 1712) 2 such that the 'dynamic' cut-off 
A is a measure for the off-shell momenta. The Hamiltonian matrix is then strictly finite: A finite 
number of sectors, and within each sector a finite number of Fock states. Now, any N x N block 
matrix equation with the Hamiltonian H can be transformed identically |2l| into a 1 x 1 block 
matrix equation with the effective Hamiltonian H c g, i.e. 



H eS \<t>) = M 2 



with \4>) = |l)(l|*/>) 



(5) 



As a consequence of working with the front form, the light-cone Hamiltonian Q and therefore also 
the effective Hamiltonian is additive in the 'kinetic energy' T and the 'effective interaction' U, i.e. 



M 



T 



U 



(6) 
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Note that the effective Hamiltonian acts only in the sector with n = 1, i.e. in the qq sector. 
The reduction from Eq.(0) to Eqs.@ or (||) is straightforward but complicated, involving only 



well-defined matrix inversions and multiplications [21]. A number of points should be emphasized: 
No unphysical Fock-space truncations are needed as they are advocated sometimes as Tamm- 
Dancoff truncations [E2L E3]. No smallness assumptions on the coupling constant is ever to be 



made. The algorithm expresses the 'higher Fock space projections' {n\ip) for n > 2 as a functional 



of (1|^>). Finally, one should emphasize that the procedure [21] is only well defined if the Fock 
space is denumerable like in DLCQ. As a result one ends up essentially with two contributions to 
the effective interaction, as displayed in Fig. p]. A more thorough presentation is presently being 
prepared. 

Every finite and well-defined Hamiltonian matrix has thus a unique and well-defined effective 
Hamiltonian matrix. In the former p5| and the present work one considers only the effective 
interaction as displayed on the left in Fig. |l| and disregards (not neglects!) the annihilation term. 
Without further arguing about that we consider the so defined effective interaction as part of our 
model. 

DLCQ was particularly useful to convert the matrix equation (^) to another matrix equation 
(||). Once this is done, one can leave the discretized case by going over to the continuum limit. In 
the continuum limit the matrix equation ^ is converted into the integral equation [25] 

V X 1 — X ' 

+ fdx'd 2 k± ,? mim2 , e , n Us^s2(x,k ± ;x',k[_) (f) s > a , (x', fcj_) . 
/j J Vx(l - x) x'(l -x') 12 

The kernel is written conveniently in terms of the interaction and the fermion currents, 
U = VC, with V(x,k ± ;x',kt) = - i 



(7) 



2vr 2 Ql ' 

and Ctllt2(x,k±;x',kj_) = u(k x , s x ) 7 M u(k[ , s[) u(k 2 , s 2 ) u(k' 2 , s' 2 ) . (8) 

The cut-off function Q' = @(x', kj_) should be a remainder, that the domain of integration is limited 
according to Eq.(Q). The four-momentum transfer Q\ is defined as the average of the quark and 
the antiquark, i.e. 

Ql = -ll = -\(k' l -k l f- l -{k> 2 -k 2 f. (9) 
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Finally, (3 is the dimensionless coupling constant. For QED it takes the value j3 = a with a ~ 1/137 

3 



For QCD with 3 colors it takes the value j3 = $a s . In either case, it depends on the momentum 



transfer due to the renormalization group. 

The kernel of Eq.(||) looks different from [^[| without being so. Here, we shall work in the 
Bj0rken-Drell convention for the spinors such that in lowest order of approximation holds j^jn ~ 1, 
as opposed to the Lepage-Brodsky convention Q were one has j^j^ ~ 4m-im2. In addition, a 



factor of 4-\/x(l — x)x'(l — x') is missing in the denominator due to a typing-error in [25]. - As 
illustrated in Fig.|], the effective interaction scatters a single quark with four-momentum k\ and 
spin-projection s\ into a quark with (A^s^), and correspondingly an antiquark from (/c2,S2) to 
(fc 2 ) s 2)- The eigenfunctions generated by Eq.(^) are identical with the probability amplitudes of 
the Fock-state expansion for 3 colors, i.e. 



W^ Al ,„M ± )M i;sl ^ 

S\,S2 



1 3 - 

with \x,k ± ;si,s 2 ) = —t= V (x, k±) (F (1 - x, -k±) |0) . (10) 

They depend on five entries: the three momenta of the quark and the two spin projections. One 
may or may not consider Eq.(|^) as the point of orign for the present model. The principal aim of 
the present work is to generate solutions or reasonable approximations thereto. 

Having been so explicit, the spin projections and their summation will be suppressed in the 
sequel as being performed implicitly. The integral equation can then be written more compactly as 

> i / r x / m ? + k\ 2 frik + k, 2 \ , , -> . f , , l9r >, 2m\m2& U(x, k±:x',k\) , . , ->. . . . 

M 2 <P(x,k ± ) = ^ + -? ±.U( X) k ± )+ dx'd 2 k' ± 1 2 v - L ' = lJ^( a / fc^) . 11 

V x 1 — x J J x(l - x) x'(l - x') 

In practice it is inconvenient to deal with the invariant mass-squared eigenvalues since binding 
effects are so much overshadowed by the large value of the total free mass M = mi + m,2. It is 
therefore suggestive to substitute the light-cone Hamiltonian H^c by an other operator H, which 
differs by the former only by a constant and a readjustment of scale, i.e. 

H LC = (mi + m 2 ) 2 + 2(mi + m 2 ) H . (12) 

The 'Hamiltonian' H will then have the dimension of a {mass) or {energy); in fact, as we shall 
see, it will have much in common with a non-relativistic Hamiltonian without being one. Since 
(mi + rri2) 2 is a Lorentz scalar, the approach continues to be frame independent. All what one does 
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is to remodel the kinetic energy by 

r = T(s, fcj.) = -; ■ r + — (mi+m 2 . (13) 

2(m 1 +m 2 ) V x 1 — x J 

Without changing the substance, the integral equation becomes then 

£^^1)=%^)^,^)+ /d 2 fcl . = r , = ^ e^^fcLs^feD^fcj). (14) 

j ^1(1 -i)r(l -rj 



The reduced mass m r is given by m r = mim^/imi + m-2), as usually. The eigenvalue of invariant 
mass-squared M 2 and the eigenvalue of the energy E are then related by M 2 = M 2 + 2ME ^ 
(M + E) 2 , distinctly different fr om non-relativistic kinematics. 

The running coupling constant does not depend on whether one calculates it perturbatively 
in the usual coordinates or by light-cone methods and behaves for sufficiently large momentum 
transfer like p7L |2| 



°-W 2 > = 33^MqW f " Q2>>K2 - <15) 



With the value of a s as measured at the Z-mass [|14j and three flavours (n/ = 3), one gets 

a s (M z ) = 0.1134 ± 0.0035 , thus k = 193 MeV . (16) 

The value of the 'QCD-scale' k is thus fixed and will not be changed during this work. But here 
is the problem: In a Hamiltonian approach like the present one, one needs the interaction not 
only at large but also at small momentum transfer, including zero. Taking the above functional at 
face value implies that hadronic interactions are repulsive for sufficiently small momentum transfer. 
This must plainly be wrong, of course, and the running must be modified. - But how? - Neither 
canonical renormalization theory nor the present DLCQ-formalism has a ready-to-use answer, yet. 
We therefore modify it in the simplest possible way, inserting a constant at the right place: 

a - (Q2) = ^ ln( Q2+ W 2 ) aSil - (17) 

Right or wrong, it does what one wants: The asymptotic behaviour is untouched and a s keeps 
its sign all the way down to Q 2 = 0. A similar step was taken by Richardson [24| for the instant 



form, but the value of k was considered a free parameter. In the sequel we shall follow Richardson 
setting a = 1, we shall however go much beyond that by including the correct fine- and hyperfine 
interaction. 
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The change from Eq.(15) to Eq.(|l7|) is the only assumption in the present work: Despite 
not completely true, one could say that 'confinement is introduced by hand' like in conventional 
potential models. The crucial difference to other models is that the correct asymptotic freedom 



1 27, £3] is combined with the virtue of having no free parameter and some further pleasant aspects: 

• The model depends only on the quark currents and is therefore manifestly gauge- invariant. 

• The model is formulated in the front form and therefore frame- and boost-invariant. 

• The usual 'recoil problem' is absent in a momentum representation. 

Finally, one should mention here already that the Light-Cone formalism has a different operator 
structure than the conventional instant form: The total angular momentum J and its eigenvalues 
cannot be used to classify the eigenstates Q. The non-initiated reader might find it quite surprizing 
that the wavefunctions can never be invariant under spatial rotations about the y- or the x-axis. 



3 Analyzing the Model by Transforming Integration Variables 

Like in other models, the flavor quark masses are free parameters subject to be determined by 
comparing with experiment. Natural candidates are the pseudoscalar (0~) and vector mesons (1~), 
which are believed to have total angular momentum L = 0. Their total spin is thus J = S. For 
two quarks, the total quark spin S can only take the value S = or S = 1. Since the flavor quark 
masses potentially range from a few MeV up to some 100 GeV, one faces a tantalizing numerical 
problem: For any given set of quark and antiquark masses, mi and m-2, respectively, one has to 
solve Eq.(0) and to perform a calculation of similar complexity as in |25[| . Actually, in addition 
to that, one is confronted with thus far unresolved, technically highly non-trivial problems due to 
a Q 4 singularity in the denominator. Moreover, of all the calculated spectrum one only needs the 
lowest eigenvalue. Only the latter will be compared with the physical meson mass M. Most of the 
calculated and thus expensive information is thus redundant, not to mention that the calculation 
will have to be repeated a zillion of times in order to fit mi and m<i to M. — Is there a simplier 
way? The following analysis will pave the way for an ultimate, possibly analytical solution. 

The fundamental equation is unpleasant to the practitioner. In Eq.(|?]) two of the three integra- 
tion variables are dimensionful (k±) and run from — oo to +oo, while the third one is dimensionless 
and has support only on the finite domain < x < 1. It is therefore more than convenient to 



9 



change variables from x to k z , such that k z has the same support and dimension as k±. This can 
be achieved by 

x = x{k z ) = ^ + + E ^ , with Ei = ^Jm 2 + k 2 + k? = ^Jmf + k 2 for i = 1, 2 . (18) 

The physical interpretation of k z is of secondary importance. One notes with pleasure that k z = 
and k± = produces the correct 'static value' x = m\j(m\ + 7712) and that one can formally define 
a three- vector k = (k z ,k±). The practitioner therefore will conclude that k z is 'something like' the 
relative momentum in z-direction. The advantages are obvious, particularly the notation is more 
compact than in the standard front form. For example, the free part of the invariant mass squared 
appearing in Eqs.(Q) and (0) becomes by direct substitution simply 

rr^+kl + rr^ + kl = + 2 
x 1 — x 

It looks as if one changes the single particle four-momentum state from the front form with k^ = 
(ki, kf, k± 1 ) and k% = (k^ , k^ , kj_ 2 ), to a state in the instant form with 

K = K k ±1 , kf) = (E 1: k) , and k% = (k° 2 , k ±2 ,kl ) = (E 2 , -k) . (20) 

The transfer and the mean of three-momentum along the quark line is then given by 

q = k-k' and p=-(k + k'), (21) 

respectively. Along the antiquark line, both have the opposite sign. But all this does not imply that 
one works in the instant form. Indeed, one has explicit residues from the front form, as follows. 
Because of boost-invariance, Eq.(^) is independent of the total longitudinal momentum P + . But 
P + still lurks in the background, since the Fock-states 'on the left' (\x, k±)) have the same P + than 
those 'on the right' (\x', kj_)). One can express this fact in instant form variables as well, using 
P + = P° + P 3 . Since P 3 = k Zj \ + £^2 = 0, one is left with P° = P ', and therefore has for every 
matrix element: 

E 1 + E 2 = E[ + E' 2 , or k 2 = k' 2 . (22) 

Because of the latter, one has also individually E\ = E[ and E 2 = E' 2 . Obviously, the interaction 
does not change the size of k, it only changes its direction: The frame-independent front-form 
formalism is equivalent to the center-of-mass frame of the instant form. This, of course, is a source 
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of great simplification. For example, the momentum transfer and the mean momentum are always 
orthogonal, 

Q-P= \{k-k'){k + k') = \{k 2 ~ k' 2 ) = , (23) 

and four-momentum transfer is always identical with the three-momentum transfer. The latter is 
verified easily. One inserts the single four-momenta from Eq.(|T9|) into the defining Eq.(||) and gets 

Ql = \ [(k[ ~ kif ~ (E[ - E,) 2 + (k' 2 - k 2 f - {E> 2 - E 2 ) 2 ] = q 2 , (24) 

This way, the interaction V depends only on the spatial momentum-transfer. 

It therefore must be our aim to express all functions appearing in the light-cone integral equation 
(|i~4l) as functions of the instant form variables k and k' , i.e. to identically rewrite it as 

Ecj)(k) = T{k)(j)(k)+ fd 3 k' @(k')A{k,k')B(k, k') C(k,k')V(k, k')(j){k') . (25) 

The notation can be chosen such that the Jacobian B and the current term C are dimensionless 
with value '1' by lowest order of approximation, as will be unrevealed gradually. Moreover, rather 
than by k and k' , we shall express them often by the above q and p, since they have a very 



special meaning in integral equations like Eq.(25), as follows. Rewrite the latter conveniently 



as Ecj)(k) = fcfik' F(q,p) (fr(k'). Multiply the whole equation with exp(ik ■ x) and subsequently 
integrate over d 3 k. Making use of the definition 

iP{x) = Jd 3 k e ilg (f>(k') , and thus F(x,p) = Jd 3 q e^ s F(q,p) , (26) 

one straightforwardly arrives at an eigenvalue equation of the Schrddinger type, i.e. 

Eip(x) =F(x,p) ip(x) , with p= -iV x . (27) 

Thus, the analogue of q is x, the position of the particle in the center-of-mass frame, and p is its 
conjugate momentum. Consider for example the interaction term in Eq. (|25|) . Since B = C = 1 will 
hold to lowest order, see below, one gets 

*«*> = -s? 2 ?^ • thus ^> = !(-; + « 2 '-)- < 28 > 

The Fourier transform generates a Coulomb potential plus a confining potential with linearly rising 
walls, to a high degree of accuracy [^4|]. This potential is plotted in Fig. || versus r = \x\. 
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Let us begin the programme with the kinetic energy T(k) = y(Ei + E 2 ) — M j /(2M). It 
is trivially a function only of k. The cut-off function 0' = 0(fe') is related to the mass-scale A 
introduced in Eq.(|j) to render the Hamiltonian matrix finite. The latter becomes now 
A 2 + (mi + m 2 ) 2 , which in turn can be cast identically into 

?2 /AA 2 A 2 + 4mxm 2 

-\2J A 2 + (mi + m 2 ) 2 " 1 j 

As expected, the cut-off is invariant under spatial rotations. Next, using the Gordon decomposition, 
the currents u(ki, Sj) 7^ u(£;-, s^) are conveniently decomposed into their convective and spin parts, 



fi ( tl , S1 )7%«,»;) = ^(1 + ^4 + ^4^) • (30) 

2mi y {L\ + rrii) z {Ei + m\) z 1 



u(ki,si)ju(k' 1 ,s' 1 ) = - — (2p — iq X <jV . , with R = iqxp. (31) 

2mi x ' 1 



Due to the validity of Eq.(22) these expressions are simpler than usually. Along the antiquark line 
one must change the sign of both k and k' and can conveniently substitute a by r in order not to 
confuse the spin of quark and antiquark. The current term C = j^j^ becomes then explicitly 



(Ei + rri!) 2 (E 2 +m 2 ) 2 J\ (E 1 + mi ) 2 (E 2 + m 2 ) 2 J 

+ 5 fMzlj^ ^-^Xf\ w . th 5 = (g 1 +m 1 )(^+m 2 ) 
V Ei + mi /V E 2 + m 2 / 4mim 2 

As one expects for a Lorentz scalar it is invariant under spatial rotations. Finally, the Jacobian 
of the transformation Eq.(|i~8|) is evaluated in terms of the dimensionless transformation function 
B(k, k 1 ), implicitly defined by 

m r 



= dx' = A(k,k')B(k,k') dk z . (33) 



y/x(l — x) x'(l — X' r 

By means of the identities 

dx _ (E l + k z )(E 2 -k z ) _ ( El + k z )(E 2 -k z ) 

dk z ~ El E 2 (E 1 + E 2 ) ^ x ^- x )- {El + E2? • ( 34 ) 

one gets straightforwardly 



Here A is the only term which explicitly 'breaks' rotational invariance since only here appears k z 
explicitly. This is no failure, but must be true for any front form Hamiltonian. As a consequence, 
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the wavefunction of the singlet state, for example, cannot be strictly invariant under rotations in 
three-space. One even has numerical evidence for such a 'violation' due to the pioneering work of 
Krautgartner pq ]: The singlet wavefunction displayed in his Fig. 11 show them clearly. First taken 
as a numerical artifact, they are meanwhile confirmed by an independent and improved calculation 
|26|| . But these 'violations' occur at a momentum scale much larger than the Bohr momentum, in 
a region where the wavefunction is down typically to 10 -3 — 10~ 4 of its peak value. They seem to 
be unimportant. 

Having defined all terms in Eq.(|i~4|), it should be emphasized strongly that its solutions must 
be identical with those obtained from the original light-cone equation (|^), of course after having 
substituted k z back by x according to Eq. (|l8|) . It is however much simpler to deal with from the 
point of view of interpretation, of approximation, and last not least from the calculational point of 
view particularly for m\ ^ m<i- 



4 The Retarded Schrodinger Equation 

The instant form has the garstly property of having square-roots scattered all over the place. 
This is particularly unpleasant when one aims at Fourier transforms. Can one find a systematic 
approximation scheme to remove them? - Once more this is possible by means of the Brodsky- 
Lapage cut-off. Indeed, one always can choose A such that 

— * < 1 , for mi < ^2 • (36) 
mf 

If this holds for the lighter particle it even more so holds for the heavier one. All square-roots can 
then be expanded safely, i.e. 

k 2 1 / k 2 \ 2 

E t ^m l + — , if « i . (37) 



2m. 



This looks like a non-relativistic approximation but the opposite is true: In the worst case, our 
choice allows for the ultra-relativistic velocities of the lighter particle up to \k\ ~ m± on the one 
hand, and on the other hand for a consistent and systematic expansion up to second order like 

A = 1 + ^ + ^ + , with a = H^l, (38) 

2m r 8m~ 2mf, \ M J m\ + 7712 

^ = 1 — 7r~~~2 _ cT~~2 ' with \ = i ft + ^) » ( 39 ) 

2m z q 8m~ m l q M \mf m^J 
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respectively. The expansion of the current term C of Eq.fl32|) and of the product of the current 
with the Jacobian run analogously and give 

B - i + gf^-t-g-^- P"/) -(^W + l-g + ^j-Li, (41) 

2 \m„ 17111712/ 4mim-2 4mj tej m\m,2 



-"^ = 1 + o n 9 J r- y- + n , 

im\m2 oniy 4mim2 4mj 4m-2 m\m,2 

respectively. For completeness we finally define S, the kinetic energy T, and 0': 

S=^(a + r), T = ^> and ©(£') = !• (43) 

The last step needs a comment. After having expanded in the manner of Eq.(|37|) one has to verify 
in principle, that the results depend on the cut-off A at most weekly. One can conjecture however 
that the wavefunction decays sufficiently fast as to act itself like a cut-off. One thus can drop A at 
this stage. All one has to do at the end of a calculation is to check the condition set in Eq.(| 
but this can and will be done by an expectation value. 



One now is prepared to Fourier transform the whole Eq.(25) into an eigenvalue equation of the 
Schrodinger type 

Hifj(x) = Eip(x) . (44) 

When transforming a momentum-space function like RV = i(q x p)V(q) one gets in a first step 

(V x p)V. By definition, the operator V can act only on the next V to the right. Since V is 

-» x dV 

spherically symmetric one sets Vl^ = — — and absorbs x into the angular momentum operator 

r or 

L = x x p. The Hamiltonian operator H in 'Schrodinger representation' (with p = —iKV) turns 



then out as 



2m r V mi +m2/ 8m 2 Amirri2 

ldVfa-L t • L S-L\ tt 
+ :^1Z2+7Z2-— T +^a. (45) 



For a spherical V holds (<r x V) • (r x VV) = §(<? • f)V 2 V and with a ■ r = 25 2 — 3 one arrives at 

m\ + m 2 (mi + i7i2) J 6m\m2 

(46) 
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The rotationally asymmetric H a is given below. We emphasize here already that the precise form 
of V need not be known in performing the step from Eq.(^) to Eq.(^). If one sticks to QCD in 
the above form, Eq.(|28|), one has 

VV(r)=/?f— + Att6(x)) and = (3 (- + \) , with (3 = % ~ 0.93 . 

\ r J r or \rr 6 J 27 

(47) 

If one sticks to QED, one chooses (3 ~ 1/137 and sets the QCD-scale n to zero. These equations 
are at the core of the present work. They are still identical with the original equations (|l4|) or 
(f25|) and include the correct retardation without smallness asumptions particularly on the coupling 
constant or the mass parameters. 

The 'Retarded Schrodinger Equation' has a wonderfully simple structure and can be interpreted 
with ease. One notes that the 'potential' V plays a different role in the different terms of the 
equation. The first term is the kinetic energy. Like in non-relatistic quantum mechanics the 
reduced mass governs the scale. Here the potential plays the role of an 'effective mass'. In the 
second term V appears in its natural role as a potential energy with an attractive Coulomb and 
a linear string potential. In the second last term, the analogue of the 'fine interactions' of atomic 
physics, the potential plays the role of a coupling constant for the spin-orbit interaction. Last not 
least, in the third and the fourth term one finds the analogues of the atomic 'hyperfine interaction' 
with the difference that they become of leading importance in the hadronic case. They contribute 
in the third term to the true potential energy and for example modify the 'Coulomb' part of V, 
while the coefficient of S 2 can be interpreted as a (color-) magnetic interaction of the two spins. 

One should note that some parts of the Retarded Schrodinger equation are a direct consequence 
of working in the front form. We have tried to account for that in the notation. In Eq.(|4l|), for 
example, one notes that the current term has no contribution in q 2 , while Eq.(f42|) does contain a 
term like that. They are caused by the explicit appearance of the light-cone momentum fraction 
x, i.e. are caused by the Jacobian BA. The asymmetry Hamiltonian H a has the same orign, and 
following Eq.(|38|) has the Fourier transform 

a z dV a 2 d 2 V dV ( 2 2m r \ p z 

H a = — - — ^-—5- + — [a 2 + . (48) 

2im r oz 8mf. az 1 oz \ M ) 2m£ 

It looks like a self-induced hadronic interaction of the quadrupole type. 

How serious is the Retarded Schrodinger Equation to be taken? - We have actually no argument 
what could have been done wrong in the above, certainly nothing on purpose or with an unjustified 
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assumption. At the very bottom the omission of zero-modes and, perhaps more important, the 
omission of the two-gluon effective annihilation graph might perhaps be influential, and has to be 
investigated in the future. But not even the 'Richardson potential' can be crucial. Small distances 
x — > correspond to large momenta in Fourier space, and at least there asymptotic freedom must 
be correct. In the opposite limit, at very large distances corresponding to small momenta in Fourier 
space, the detailed form of the running coupling constant might possibly matter. But there, the 
bound-state wavefunction is small from the outset and possible effects are unlikely to become 
important. Not even the consistent expansion to second order in the momenta is an assumption, 
since it can and must be checked a posteriori by 



One thus has to live with these equations. Rather than to agonize further on these almost philo- 
sophical matters, to which we shall come back in the summary, we proceed now by asking for 
practical consequences. 

5 Solving the Model by Parametric Variation 

It will take some time and effort to work out all the many consequences of the Retarded Schrodinger 
Equation (f46|), which as mentioned should hold both for QED and QCD. In the sequel, we therefore 
shall restrict ourselves to calculate only the ground-state masses of the pseudoscalar and vector 
mesons. If one leaves aside the recently discovered top-quark t and restricts to 5 flavors (u,d;s,c;b), 
one has thus 30 different physical mesons, since charge-conjugate hadrons have of course the same 
mass. 

One cannot calculate these masses, however, without knowing the mass parameters mi and 
rri2 of the quarks. These cannot be measured in a model-independent way, although some people 
pretend one can do so. Here, we shall adopt the point of view that the quark masses have to 
be determined consistently within each model, for the better or the worse. One has thus 5 mass 
parameters to account for 30 physical masses. Which ones should be selected to fit? — There are 
142 506 different possibilities to select 5 members from a set of 30, and we have to make a choice: 
We choose the five pure qq-paiis. Even that is not unique: Shall one take the pseudoscalar (0 - ) or 
the vector mesons (1~)? - We shall take both, and compare the results. 




(49) 
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But doing so, one runs into the problem of choosing the chiral composition of the physical 
hadrons. One would be rather surprized if the pions, for example, could be understood as eigenstates 
in a simple potential model like Retarded Schrodinger Equation (|4(f). Thus, in the aim to avoid 
hell's kitchen by choosing a certain flavor composition, we shall do even worse and replace the 
'pions', for example, by 'quasi-pions' with the same physical mass: All our hadrons shall be pure 
gg-pairs - by fiat. For example, we shall identify the uu-, ud-, du-, or dd-eigenstates with the 
quasi-7r°, quasi- 7r + , quasi- tt~ , or the quasi-77, and so on. These 'crimes' can be revoked easily in 
future though not the present work. They are by no means compulsory to the model, as we shall 
see. 

Our problems lie in another ball park. It is clear that we cannot deal head-on with the full 
complexity of the Retarded Schrodinger Equation (fl6|). Which part of the Hamiltonian should 
one select in a first assault? - Some help is gained by the rather unique property of the light-cone 
Hamiltonian: Kinetic energy and interaction are additive, and so is the Hamiltonian. Because of 
this additivity one always can select an 'interesting part' Hq, i.e. H = Hq + AH and check a 
posteriori, by calculating ( AH ) in perturbation theory, whether the selection makes sense, or not. 

Since the scalar and pseudoscalar mesons have probably no orbital excitations, i.e. are primarily 
s-waves, one can disregard the spin-orbit couplings and choose 

B ,„( l+ J^)£ +vlT) -p'(-^-<™-™r) + *^8>. ( 50) 

V mi +7712/ 2m r amf. 1 mi + m2 \jn\-\-m2) I 3mim2 

Even that looks to complicated for a start-up. We therefore select those terms which in one way or 
the other have turned out to be important in past phenomenological work, namely a simple central 
potential plus the triplet-interaction mediated by the total spin. Our first choice is therefore 

rr P 2 ~kt I \ 2K & P 2 P 2 * 2 & a 2 

H = ^— + V(r) + = £- - tL + + I3n 2 r . (51) 

It is always possible to choose a spinorial representation in which S z and S 2 become diagonal. In 
this case, the latter can be replaced by the eigenvalue S e = S(S + 1), taking the values and 2 for 
the singlet and the triplet, respectively. 

How does the wavefunction for the lowest state look like? - For a pure Coulomb potential, i.e. 
for k = 0, the solution has the form 

il>(x) = -^Aie" Ar . (52) 

\/7T 
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Table 1: The Eavor quark masses in MeV, as obtained from a fit to Eq.(57). — The first row refers 
to a fit for the singlets, the second to the one for the triplet. 



Flavor mass 


u 


d 


s 


c 


b 


From fit to 0~ 


2.3 


155.6 


430.6 


1642.3 


5330.8 


From fit to 1~ 


222.8 


236.2 


427.2 


1701.3 


5328.2 



One the other hand, ommitting the Coulomb part, a linear potential be solved in terms of Airy 
functions or its integral transforms j2j|. If one has both, one will have some mixture of the two. 
In a numerical calculation, these can be found with almost arbitrary accuracy. But for the present 
start-up check even that requires too much effort. 

Rather we shall pursue a variational approach. We choose Eq.(|52"|) as a one-dimensional param- 
eter family, subject to determine the one free parameter A. Of course, one could take other families 
as well, like for example harmonic oscillator eigenstates, but with the above all expectation values 
are particularly simple, i.e. 

(V| f | V) = A 2 , (VI \ IV) = A , and (V| r |V) = ^ • (53) 

That's all what one needs to calculate the expectation value of the energy 

— i 1 1 r i / \ A 2 2S e K 2 3/3k 2 1 , s 

Since we deal only with a variational approach to the ground states, we are not in conflict with 
the statement that the wavefunction cannot be purely Coulombic. For the pure Coulomb case, 
the 25- and IP-states would be degenerate and the ratio |V2s(0)| 2 /|Vis(0)| 2 = 0.125, as opposed 
to the values ~ 0.63 and ~ 0.50 extracted from experiments on charmonium and bottomium |i~5]| , 
respectively. 

Since we address ourselves to calculate the total invariant mass of the hadrons we return to the 
front form invariant mass-squared operator Hlc = M 2 , i.e. to M 2 = (mi + mj) 2 + 2(mi + rri2)E. 
Specializing now to equal masses m\ = mi = m one preferably converts the variational equation 
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Table 2: The validity check. — The first row displays the values of (X/m) 2 as obtained from 
the mass fit to the singlets, the second row those from the mass fit to the triplets. — Note: If 
(X/m) 4 > 8, the solution has to be rejected as a consequence of Eq.filQ), since (p 2 ) = X 2 . The 
extremely large value for the u-quark in pseudoscalar fit gives a good example for such a case. 



A 2 
m 2 


u 


d 


s 


c 


b 


For fit to 0" 


293 


1.45 


0.53 


0.25 


0.22 


For fit to 1" 


0.79 


0.76 


0.48 


0.25 


0.22 



([54]) in terms of the dimensionless variables 

A ^ m fM\ 2 , . 

s = -, £ = - , and W = — . 55 

We measure thus all energies and masses in units of the one fixed scale in our problem, the QCD 
scale k. The variational eigenvalue equation reduces then to the handy expression 

At fixed values of the parameters (£, /3, S e ), the variable A and thus s must be varied such that the 
energy (or the mass) become stationary, i.e. 



dW(s;0 



Os 



M 



2k 



2 



0, thus — =W(s*(£))=W*(0. (57) 



=**(£) 

This leads to a cubic equation in s which can be solved analytically in terms of Cardan's formula. 
In special cases they can well be approximated by a quadratic equation, namely when s* > 1 or 
when s*<l. We got accustomed to refer to these two regimes as the Bohr and the string regime, 
respectively. In the Bohr regime the Coulomb potential dominates the solution and the linear string 
potential provides a correction. In the string regime the linear string potential dominates, with 
the Coulomb potential giving a correction. Solutions in the string regime, however, imply that the 
ratio (p 2 ) /m 2 = (X/m) 2 becomes so large that one is in conflict with the validity condition Eq.(fi"9"|). 

Rather than to give explicitly the straightforward but cumbersome formalism, we present the 
analytical results in the graphical form of Figure |3|, both for the singlet and the triplet. As displayed 
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there, the total mass M = 2k\/W* is almost linear in the quark mass, with small but significant 
deviations. In line with expectation, the hyperfine splitting decreases with increasing quark mass. 
Less expected was that the splitting increases so strongly with decreasing quark mass. As one 
can see in the Figure, the mass of the singlet states takes off from the value at vanishing mass 
almost like a square-root as opposed to the triplets which start off at a finite value. Therefore, if 
one determines the mass of the li-quark by fitting to the quasi-pion mass, as displayed in Table |], 
the value of s becomes extremely small. Consequently, one is in the string regime and the u-quark 
becomes ultra-relativistic. The validity condition is then badly violated, as compiled in Table ||. 
It is worth noting that the value of the up-quark mass is rather close to what is reported for the 
so called 'current mass'. For the r\ and the rf, the scaling variable s is of order unity, while for 
the quasi-r/c one definitely is in the Bohr regime. Here the masses are similar or close to what is 
referred as the 'constituent-quark' mass. 

Since singlet and triplet are so close for s>l and since nobody wants to fit the pion anyway, 
quasi or not, one fits the quark masses preferentially with the vector mesons. In the lack of 
empirical data we have set M„ b = Mx, which should be of minor importance in the present 
model, see Figure |3[ As seen from Table |l|, the flavor masses are now in close agreement with 
the constituent quark masses, particularly the up- and down-quarks. As compiled in Table ||, the 
smallness condition is now satisfied excellently, everywhere. 

Obviously, it would be a minor modification to introduce some preferred chiral structure of the 
light mesons, and to drop the simple and simple-minded construction of quasi-mesons. But we have 
an other objection, here. We ask: With the so fixed flavor quark masses, how bad or how good are 
the remaining 25 meson masses described? - The procedure runs quite analogously, except that it 
is now easier. With fixed masses, one varies A in Eq.(^) such that the energy becomes stationary. 
Of course, this has to be done separately for each flavor composition. We compile the results in 
Table |||. The table contains also a comparison with the experimental values, to the extent the 
latter are known. It is clear why one should exempt the quasi-pions from a detailed comparison 
with data, but for the rest it was not really anticipated that they would agree on the level of about 
one part in hundred. Some of these mesons have not yet been found. The present model predicts 
then for example 

M{Bf) = 6494 MeV , M(B* ± ) = 6501 MeV . (58) 
Such a table was not given before. In all due respect for the work with lattice gauge theory, the 
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agreement of the above equations with data is not so much worse, particularly in view of their 
simplicity and the absence of any free parameters. 



6 Summary and Conclusions 

Applying DLCQ to QCD in the light-cone gauge A + = and disregarding zero modes, one derives 
first an effective interaction acting in the qq space. We emphasize that no assumption on a possible 
smallness of the coupling constant has to be made, nor that it is necessary to truncate the Fock 



space in the manner of Tamm [|22|] and of Dancoff [23]. Gauge invariance, therefore, is not violated. 



Once this is achieved, one can drop discretization and convert the matrix into an integral equation, 



Eq.(ll). Unfortunately we cannot be more explicit in the present context. We therefore consider 



this equation as part of the model and the point of departure in the present work. - We discover 



next a simple substitution of integration variables, Eq.(18), which allows to transmute the equation 
identically into Eq.(^), in which appear only the conventional instant form variables like the usual 
momenta and energies. One wishes to Fourier transform the still front form equations to usual 
space coordinates but is hampered by the 'square roots'. A poor man's solution is their consistent 
expansion up to second order. The prize to pay is Eq. (|49|) , i.e. a smallness requirement a posteriori 
which has to be met in the solution. One ends up with a differential equation in configuration space, 
Eq.(|50|)) which looks like a non-relativistic Schrodinger equation without being one. It contains 
the full retardation of a relativistically correct front form equation, and therefore is referred to as 
the Retarded Schrodinger Equation. Its solutions can be Fourier transformed back to momentum 



space and subjected to the transmutation equation (18), to get the wavefunction expressed in light- 
cone momenta particularly longitudinal momentum fractions, which are needed so urgently for the 
theoretical predictions of experimental high-energy cross sections, and of structure functions and 
the like. 

One should emphasize that the Retarded Schrodinger Equation has no free parameter: The 
coupling constant and the quark masses have to be determined from the experiment. Like in the 
familiar quantum electro-dynamics one has no way to calculate them. Fitting the strong coupling 
constant at the Z-mass and the 5 quark flavor masses to 5 selected vector mesons fixes all our 
freedom. The rest is structure: The 20 remaining pseudoscalar or vector masses are then predicted 
and presented in Table 0. In comparison with the experimental data they are not much worse 



than those from conventional phenomenological models [15, Pul, or from heavy quark symmetry 
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1 17], or even from lattice gauge calculations Jig, I3 20]. The pions, of course, remain mysterious 
particles like in every other model not specially designed for them. One wishes now to solve the 
full Retarded Schrodinger Equation ( |50|) rather than only its caricature in Eq,([5l|). 

Conclusion: If such a poor model can do so well one must be on the right track. It seems that 
the front form and DLCQ have made a big step forward. Unfortunately, it took ten years but one 
has been faced with all the problems of formulating a Hamiltonian approach to gauge field theory, 
in addition to the problems inherent to the front form. 
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Figure 1: The effective interaction in the qq sector. — By absorbing (or emit- 
ting) an 'effective gluon', a single quark state with four- momentum k\ and spin- 
projection s\ is scattered into the quark state (k^s'i)- Correspondingly, the 
antiquark is scattered from (&2,S2) to (k^s^). 
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Figure 2: The quark-antiquark potential V(r). Inserted is also the Coulomb 
(dashed-pointed) and the confining potential (pointed). 
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Figure 3: Bound states of qq-pairs versus quark masses. — All masses are given in units of the 
QCD-scale k. The upper curve refers to the triplet (S e = 2), the lower to the singlet (S e = 0). The 
masses of some vectormesons (p°,u>,(f),J/ijj) are marked by (A), those of the some pseudoscalars 
(k°,V,v',Vc) by (*). 
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Table 3: The masses of qq-hadrons are compared with experimental values. - Flavor quark masses 
used are inserted in column 2. They in turn come from a fit to the vector mesons, as discussed in 
the text. - Within each box, the first line refers to the hadronic symbol of the meson; the second 
line gives the calculated (measured) vector mass in MeV; the third line accounts for the calculated 
(measured) pseudoscalar mass in MeV, and finally the fourth line specifies the hadronic symbol of 
the pseudoscalar meson. - Note: The smallness parameter (A/to) 2 satisfies always the condition 
set in the text. 





m q 
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c 


b 


u 


222.8 


P° 

768(768) 
714(135) 


P + 
773(768) 


K*+ 
910(892) 


2110(2007) 


B *+ 
5712(5325) 


d 


236.2 


658(140) 


to 

668(549) 

7] 




D*~ 
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427.2 


825(494) 
K~ 


831(498) 
K° 


1019(1019) 

953(958) 
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D*- 
2156(2110) 


Bf> 
5735( — ) 
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1701.3 


2079(1865) 
D° 


2078(1869) 
D+ 


2131(1969) 
Df 


3097(3097) 
3082(2980) 

Vc 


B* + 
6502( — ) 
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5328.2 


5701(5278) 
B~ 


5698(5279) 
B° 


5726(5375) 


6495( — ) 
Be 


T 

9460(9460) 

9455( — ) 
Vb 
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